Emergent quantum confinement at topological insulator surfaces 
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Bismuth-chalchogenides are model examples of 
three-dimensional topological insulators. 1 Their 
ideal bulk-truncated surface hosts a single spin- 
helical surface stateJSES which is the simplest 
possible surface electronic structure allowed by 
their non-trivial Z2 topology. 4 They are therefore 
widely regarded ideal templates to realize the pre- 
dicted exotic phenomena and applications of this 
topological surface state. However, real surfaces 
of such compounds, even if kept in ultra-high vac- 
uum, rapidly develop a much more complex elec- 
tronic structure 5 6 whose origin and properties 
have proved controversialM*"2JII Here, we demon- 
strate that a conceptually simple model, imple- 
menting a semiconductor-like band bending in 
a parameter-free tight-binding supercell calcula- 
tion, can quantitatively explain the entire mea- 
sured hierarchy of electronic states. In combi- 
nation with circular dichroism in angle-resolved 
photoemission (ARPES) experiments, we further 
uncover a rich three-dimensional spin texture of 
this surface electronic system, resulting from the 
non-trivial topology of the bulk band structure. 
Moreover, our study reveals how the full surface- 
bulk connectivity in topological insulators is mod- 
ified by quantum confinement. 

Topological insulators (TIs) are an exotic state of 
quantum matter, guaranteed to have metallic edge or 
surface states due to an inverted ordering of their bulk 
electronic bands. The corresponding topological invari- 
ants dictate that there must be an odd number of such 
states intersecting the Fermi level between each pair of 
surface time-reversal invariant momenta (TRIM). In the 
most widely investigated bismuth-chalchogenide family 
of TIs, there is just one of these so-called topological sur- 
face states (TSSs) cr eating a si ngle Dirac cone around the 
Brillouin zone centreP- 3 * 5 * 19 * 20 * However, exposure to even 
the minute amount of residual gas in an ultra-high vac- 
uum chamber induces a drastic reconstruction of the sur- 
face electronic structure (Fig. [T]). The Dirac point shifts 
to higher binding energies, indicating more electron-rich 
surfaces! 5 * 6 * 8 * More importantly, additional pairs of two- 
dimensional, almost parabolic states emerge in the vicin- 



ity of the bulk conduction bandP which develop large 
Rashba-type splittings^ while new ladders of M-shaped 
states are created in the original bulk valence bandsP^ 
Given their propensity for formation, these must be con- 
sidered part of the intrinsic electronic structure of the 
surface of any realistic TI, for example when exposed to 
air or interfaced to another material. However, their ex- 
istence is not predicted by idealized theoretical models 
of the bulk-truncated surfaced and has proved contro- 
versial. Recent proposals attribute their formation to 
a variety of electronic or structural modifications of the 
crystal hostP^ but to date no single model has been 
able to simultaneously reproduce the number and binding 
energies of the experimentally-observed electronic states 
as well as the magnitude of their measured Rashba-type 
splittings. 

In the following, we show that this can be achieved 
taking account of only a single electronic phenomenon 
- a near-surface electrostatic potential variation. Such 
intrinsic electric fields are familiar from surfa ce sp ace- 
charge regions in conventional semiconductors 1211221 anc j 
are also formed at interfaces where they can be modu- 
lated using applied gate voltages. Thus, our findings are 
not only relevant to understand the free surfaces of topo- 
logical insulators, but also their use in transistor-style 
or other thin-film device architectures, which may ulti- 
mately provide a route towards room temperature spin- 
tronics. Here, we consider Bi2Se3 as a model example, 
but our computational approach can be freely adapted to 
virtually all TIs, such as other emerging Bi-chalcogenides 
or the half-Heusler alloysp^ where we would expect a 
similarly striking role of band bending at free surfaces, 
interfaces, or in devices. 

We implement such a band bending scenario within 
a 120 quintuple-layer (QL) tight-binding supercell cal- 
culation (see Methods). Without incorporating any 
band bending, our supercell model yields a single TSS 
spanning from the bulk valence to conduction bands 
(Supplementary Fig. SI), in good agreement with both 
first-principles slab calculations^ and our experimental 
measurements of the pristine surface (Supplementary 
Fig. S2). Additionally including an electrostatic poten- 
tial variation, we find a hierarchy of electronic states 
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FIG. 1: Band-bending-induced surface electronic 
structure of Bi2Se3. (a) A near-surface electrostatic poten- 
tial variation forms a surface quantum well. This shifts the 
Dirac point of the topological surface state (TSS) further be- 
low the Fermi level and causes the conduction band states to 
be restructured into ladders of Rashba-split two-dimensional 
subbands, with envelope wave functions (dark red) peaked 
close to the surface. The finite valence band width provides 
simultaneous confinement of the valence band states, caus- 
ing these to also become quantized into ladders of M-shaped 
states. These features can all be seen in ARPES measure- 
ments of the surface electronic structure (b) and are well 
described by the surface projection (i.e. real-space projec- 
tion onto the first 3 QL) of an ab-initio-denved tight-binding 
model (c) of the ideal Bi2Se3 bulk structure subject only to a 
perturbing electrostatic potential close to the surface. 



that is in excellent agreement with our experimental 
angle-resolved photoemission (ARPES) measurements of 
a doped TI surface (Fig. [ljb,c)) This agreement is re- 
markable, given that we do not adjust any parameters in 
our calculation. All tight-binding matrix elements are ob- 
tained from an ab-initio calculation of the bulk electronic 
structure. The magnitude of the potential change is fixed 
by the experimentally-measured shift of the Dirac point 
between a pristine cleaved surface and the n-type doped 
one (see Supplementary Fig. S2), and its functional form 
follows from a solution of Poisson's equation.^ We stress 



that our approach does not require any additional modi- 
fications beyond an electrostatic potential , such as a huge 
increased van der Waals gap siz d 12 * 13 * 15 * 17 *, other artificial 
detachments of a finite number of crystal layers^, or the 
formation of new topological surface states^. 

Rather, all of the new electronic states observed ex- 
perimentally arise simply due to the perturbative ef- 
fect of the electrostatic potential in the system, which 
causes a triangular-like downward bending of the bulk 
electronic bands close to the surface by almost 250 meV 
(Fig. [2](b)) . Combined with the large potential step 
at the material/vacuum interface, this creates a near- 
surface quantum well for electrons. The bulk conduction 
band states are consequ ently re structured into ladders of 
multiple subband statesP^filSZI Q ur layer-resolved calcu- 
lations (Fig. |2| show that the lowest such subband, lying 
deep within the potential well, is localized in the topmost 
few quintuple layers (QL) of the structure: The electrons 
that populate this subband are free to move parallel to 
the surface, but are strongly confined perpendicular to 
it, forming a two-dimensional electron gas (2DEG). This 
state exhibits a strong Rashba spin-orbit splitting, whose 
calculated magnitude reproduces our experimental obser- 
vations. Higher-lying states of the subband ladder are 
less strongly bound within the surface quantum well and 
contribute significant weight as far as 10-15 QL below 
the sample surface. The smaller spin splitting observed 
experimentally for these higher-lying states is again re- 
produced by our model calculations, confirming that the 
Rashba splitting is driven by the potential gradient of the 
confining electrostatic potentialP Similar to the conduc- 
tion band 2DEG, a ladder of subbands is also observed 
in the valence band near the zone center. These hole- 
like states become quantum confined between the surf ace 
and the upper edge of a projected bulk band gari^, as 
shown schematically in Fig. 1(a). We also find a small 
spin splitting of these states, which is not resolved exper- 
imentally. Its location, around the top of their M-shaped 
dispersions, indicates a subtle interplay of spin-orbit in- 
teractions between the valence and conduction band sub- 
bands.^ Thus, despite its conceptual simplicity, the band 
bending has dramatic consequences for the surface elec- 
tronic structure. Driven by an emergent role of quantum 
size effects, it leads to a hierarchy of electronic states in 
both energy and spatial extent, spanning from the two- 
dimensional limit of the TSS and lowest quantum well 
states of the conduction and valence bands to the three- 
dimensional continuum. 

Intriguingly, we find an anti-correlation between this 
electronic dimensionality, and that of the states' spin tex- 
ture. This is revealed experimentally via an asymmetry 
in the matrix element for photoemission when the sys- 
tem is excited by right and left circularly polarized light, 
respectively. Such circular dichroism has previously been 
suggested to directhJ22l or indirectly^ probe the spin of 
the isolated TSSpSEH Here, we adopt this concept to 
show that the full three-dimensional spin texture of the 
entire low-energy surface electronic structure can be ex- 
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FIG. 2: Confinement of the electronic states, (a) and (b) show the depth-dependence of the carrier density and potential 
profile near the surface of Bi2Se3, respectively, (c) Layer-projected supercell calculations reveal the varying spatial confinement 
of the different electronic states. In particular, the topological as well as the lower 2DEG and M-shaped quantum well states 
are all strongly confined within the first 3-4 quintuple layers (QLs) below the surface. 
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FIG. 3: Three-dimensional spin texture of the topological and 2DEG states. (a,b) Circular dichroism in constant 
energy contours and E vs. k dispersions measured by ARPES encode the spin texture of the topological and 2DEG states. CD 
along Y-M (b), containing only the in-plane component of the spin texture, revealing opposite helicity for consecutive Fermi 
surface sheets, (c) Quantitative analysis of the angular (0) dependence of the dichroism around the Fermi surface, shown for 
the outer Eo 2DEG state. A large sin(3<^>) contribution, also observed for the TSS (see Supplementary Information), reveals 
a significant out-of-plane spin canting of the larger Fermi surface sheets, consistent with the calculations shown in (d). In 
contrast, all other states largely retain the in-plane spin texture characteristic of classic Rashba systems, all the way up to 
the Fermi energy. Together, this leads to a rich three-dimensional spin-texture of the surface electronic structure of TIs, as 
summarized for the TSS and lowest (Eo) Rashba-split subband of the 2DEG in (e). 



4 




FIG. 4: Connectivity of the surface electronic spectrum, (a) Full electronic structure calculation along the T-M 
direction. At the F point, Kramer's degeneracy requires that the FSS is spin-degenerate (the Dirac point, DP) and each 2DEG 
state is degenerate with its spin-split pair. At the M point, this degeneracy occurs between the topological state and outer 
branch of the 2DEG states, and then between each successive pairs of 2DEG states, as shown schematically in (b). Between 
these two time-reversal invariant momenta (FRIM), no fundamental symmetry dictates degeneracy of the states, and indeed 
these all remain non-degenerate up until at least the Fermi level, visible in high-resolution ARPES measurements (c). 



tracted in a similar way. Along the T-M mirror direction, 
symmetry requires that the spin must lie entirely within 
the surface planeP^ Our measurements of CD extracted 
along this direction (Fig. |3jb)) are therefore representa- 
tive of the in-plane spin texture, and indicate an alternat- 
ing left-right-left-right-left helicity for consecutive Fermi 
surface sheets of the TSS and first (Eo) and second (Ei) 
quantized conduction subbands. Away from this direc- 
tion, dichroism of the inner Eo and both Ei 2DEG Fermi 
surface sheets vary sinusoidally with angle, (f>, around 
the Fermi surface (shown in Supplementary Information) . 
This reveals a model Rashba-like 3? helical spin texture 
of these states, with the spin polarised entirely in the 
surface plane at all energies. However, for the more 
spatially-localized TSS and outermost 2DEG states, the 
dichroism additionally develops a strong sin(3</>) compo- 
nent which is maximal along T-K and zero along T-M 
(Fig. (3^c) and Supplementary Information). This re- 
flects a strong out-of-plane spin polarization close to the 
Fermi level, which is correlated with hexagonal warping 
of the electronic statespS and again switches sign be- 
tween neighbouring states. Thus, not only does the band 
bending promote the creation of additional Fermi sur- 
face sheets, it also leads to complex band- and binding- 
energy dependent three-dimensional spin-textures of the 
resulting electronic system, as summarized in Fig. [3|e). 
Together, this will cause significant complications for the 
interpretation of spin-dependent transport experiments, 
and inter-band processes will weaken the topological pro- 
tection against backscattering commonly assumed for the 
TSS of these compounds. Nonetheless, we note that our 
resolution limited linewidths of < 0.006 A" 1 (< 0.13°) 
indicate that the intrinsic scattering remains relatively 



weak for all Fermi surface sheets. 

The non-trivial spin texture of these states also raises 
questions over their interplay. For any system with TR 
symmetry, the electronic states must be spin-degenerate 
at TRIM, even though a breaking of inversion symmetry 
by the confining potential and surface allows a lifting of 
the spin degeneracy at arbitrary fc-points. The hexago- 
nal surface Brillouin zone of topological insulators such 
as Bi2Se3 contains TRIM fc-points at the zone-centre (r 
point) and at the side-centre (M) point. Fig.|4^a) shows 
the electronic structure calculated along the entire T- 
M direction. At the T point, each Rashba-split pair of 
2DEG states clearly become degenerate. For the TSS, 
this initially appears to be complicated by the presence 
of bulk states, with the Dirac point looking to be buried 
deep in the bulk valence band. However, locally, the TSS 
is still situated within a projected band gap of the (bent) 
bulk bands, and so remains a well defined surface state 
(unlike the 2DEG states which are located within, and 
indeed derived from, the projected bulk-like bands). In- 
deed, our layer-projected calculations (Fig. [2]) reveal that 
the Kramer's degeneracy of the TSS is preserved in the 
presence of the band bending potential, contrary to pre- 
vious suggestions of a band gap opening in the surface 
spectrum at the Dirac point Q3D3I 

At the M-point, our calculations reveal that the TSS 
connects to the outer branch of the lower 2DEG state. 
Assuming a spin-down state for the TSS, the pairing 
2DEG state has no other choice but to be spin-up. Its 
spin-down partner must then become connected with 
the outer branch of the second subband to satisfy TR 
symmetry for these states, thereby dictating the overall 
up-down-up-down-up spin configuration discussed above. 
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Fundamentally, this pair-exchange follow from the non- 
trivial Z2 topology of the bulk band structure. This re- 
quires that the TSS must connect from the bulk valence 
to conduction bands. Now, however, the original bulk 
bands have become quantized into two-dimensional sub- 
band states in the vicinity of the surface, and our cal- 
culations reveal that it is to these subbands which the 
topological state connects. We note that this requires 
the outer and inner branches of the 2DEG states to di- 
verge from each other at large fc-vectors, in stark contrast 
to conventional Rashba systems^ in which the spin-split 
bands merge again at zone-edge TRIM. 

Hexagonal warping of the electronic states can ad- 
ditionally indu ce ac cidental degeneracies along high- 
symmetry lines J2SJ3B] as seen m pjg rjj However, as these 
are not protected by any fundamental symmetry, intro- 
duction of perturbations in the system can shift, or even 
completely remove, such accidental degeneracies. Exper- 
imentally, the inner branch of the lower 2DEG state and 
outer branch of the upper 2DEG state can be seen to 
approach each other at low binding energies (Fig. Qc)), 
consistent with their proximity to the accidental cross- 
ing predicted in our calculations. However, contrary to 
previous reports^lA^I, our high-resolution measurements 
reveal that the occupied bands remain nondegenerate and 
exhibit rather parallel dispersion close to the Fermi level, 
which could be indicative of an avoided crossing. Under- 
standing whether interactions not included in our theo- 
retical treatment play a role in avoiding the accidental 
degeneracies will be crucial to unravel surface-bulk cou- 
pling in these materials, and the role of interactions in 
limiting mobilities of spin-polarized surface transport in 
topological insulators. 

Methods 

Theory: Relativistic electronic structure calculations were 
carried out within the context of density functional theory 
using the Perdew-Burke-Ernzerhof exchange-correlation func- 



tional and the augmented plane wave plus atomic orbitals 
(APW-LO) method as implemented in the WIEN2K pro- 
gram.^ For all the atoms, the muffin-tin radius Rmt was 
set to 2.5 Bohr and the maximum modulus of reciprocal 
vectors K max was chosen such that RMTK max = 7.0. The 
primitive cell was considered to be hexagonal (space group: 
R3m) with lattice parameters and atomic positions taken 
from experiment.^ The corresponding BZ was sampled using 
a 10 x 10 x 3 fe-mesh. To simulate the effect of band bending, a 
120 QL tight-binding supercell Hamiltonian was constructed 
by downfolding the APW-LO Hamiltonian using maximally 
localized Wannier functions (MLWF's) P^"^ We chose valence 
p orbitals of Bi and Se as the projection centers of MLWF's. 
The bending potential was obtained by solving the coupled 
Poisson-Schrodinger equation,^ assuming a static dielectric 
constant of 70 P The resulting potential was then added to 
the onsite terms of the tight-binding sucpercell Hamiltonain. 
The (QL-projected) surface band structure was eventually ob- 
tained by diagonalizing the supcercell Hamiltonian (and pro- 
jecting the resulting eigenstates onto the MLWF's of respec- 
tive surface QL's). 

Experiment: Bi2Se3 crystals were prepared by reacting high 
purity elements (5N) of Bi and Se in sealed, evacuated quartz 
tubes at 850°C for 2 days, followed by cooling at 2-3°C/h to 
650°C. The crystals were annealed at this temperature for 7 
days before quenching to room temperature. ARPES mea- 
surements were performed at 10 K using circularly-polarized 
light with photon energies of 20-30 eV and a Scienta R4000 
hemispherical analyser at the SIS beamline of the Swiss Light 
Source. Surface doping was achieved by depositing potassium 
on freshly-cleaved samples from a properly outgassed SAES 
K getter source. 
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SI 1. ELECTRON DOPING OF THE SURFACE 

The surface projection of our ab-initio-derrved tight- 
binding model for the ideal bulk-truncated surface of 
Bi2Se3 is shown in Fig. |Sl[a). Besides the projected 
bulk bands, the topological surface state gives rise to a 
clear Dirac-like feature spanning from the bulk valence 
to conduction bands, in good agreement with previous 
idealized first-principles models of the surface electronic 
structure of topological insulators!^ A very similar elec- 
tronic structure is also seen in our experimental mea- 
surements of the pristine surface of Bi2Se3, performed 
immediately after cleaving the sample in ultra-high vac- 
uum. Additionally imposing a near-surface electrostatic 
potential in our model (i.e., a downward band bending 
resulting from surface doping) leads to a drastic recon- 
struction of the electronic structure, as discussed in the 



main text and shown for comparison in Fig. SI [b). 

Surface doping via adsorption of residual gas from the 
vacuum onto the surface of Bi2Se3, and related com- 
pounds, is already sufficient to drive time-dependent 
changes of the electronic structure, which ultimately 
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FIG. SI: Projected band structure from a 120QL ab-initio- 
derived tight- binding model of the ideal Bi2Se3 structure (a) 
without and (b) including an additional near-surface band 
bending. 
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FIG. S2: ARPES spectra taken from the pristine (left) and 
K-dosed (right) surface. The Dirac point shifts downwards 
by 240 meV, accompanied by the quantization of the bulk 
conduction and valence band states. 



leads to such surface quantum well format ion 1213] jj ow _ 
ever, creation of such band bending can be drastically ac- 
celerated by deliberately adsorbing minute quantities of 
a range of molecules or atoms on the surface Here, we 
employ alkali metal dosing, depositing <C 0.1 ML (from 
XPS analysis) of potassium on the freshly-cleaved sample 
surface from a properly outgassed SAES K getter source. 
This causes the spectral changes shown in Fig. |S2| includ- 
ing a shift of the Dirac point to higher binding energies by 
240 meV, and the emergence of Rashba-split conduction 
band subbands and M-shaped valence band subbands, 
as also seen in our theoretical model (Fig. |Sl~|b)) and 
discussed in detail in the main text. 



SI 2. QUANTITATIVE ANALYSIS OF 
CIRCULAR DICHROISM 



The experimental geometry for the circular dichroism 
measurements is shown in Fig. S3 'a) . The analyser en- 



try parameter 



trance slit is oriented along the T-M direction, and the 
full electronic structure is measured by tilting the sample 
around this axis. A clear intensity difference can be ob- 
served between the Fermi surface measured in this way 
using right and left hand circularly polarized light, re- 
spectively (Fig. [S3|^b-d)). Such circular dichroism (CD) 
has recently been shown to encode the spin-texture of 
the topological surface state (TSS) in topological insula- 
tors (Tls)P^We use this method to additionally extract 
the spin texture of the 2DEG states, as sumarized in the 
main manuscript. We note that, for the relative changes 
in CD between the bands that we study here we can ne- 
glect final-state effects, as transitions from all of these 
closely-spaced initial states will go into the same (broad) 
final state. 

To enable quantitative analysis, we define an asymme- 



7 c +(^)-J c ~(0) 



where I c± ((f)) is the normalized intensity of photoemis- 
sion measured using right or left circularly polarized 
light, respectively, at an azimuthal angle <f) around a 
Fermi surface contour (defined in Fig. S3 |e)), with a 
node in the CD always found for (f> — 0. Such asym- 
metry measurements are shown at the Fermi level in 
Fig. S3 f) and at binding energies of 0.15 eV and 0.4 eV 
in Fig. S4[c) and (f), respectively. These are fit by the 
function Asm(<fr) + Bsin(30), with the A, B coefficients 
shown in the figures. CD varying with a simple sin(0) 
dependence around the Fermi surface arises from the in- 
plane winding of the spin around the Fermi surface, while 
the sin(3</)) component r esults f rom coupling to the out- 
of-plane spin component EEJIL21 The ratio of the A and 
B components in the fit therefore gives a good measure 
of the degree of out-of-plane spin polarization. We note 
that this is the minimal model consistent with the C3,, 
symmetry of the surface^. Negligible improvement in 
the fits was found when including possible higher-order 
terms (such as a sin(60) componenli 1 1 1 13 l) , so we neglect 
these in the following. 

The asymmetries around the Fermi surface are shown 
in Fig. S3 i). The asymmetry changes sign upon moving 
inwards from the outermost Fermi surface sheet of the 
topological state to the subsequent Fermi surfaces of the 
lowest Rashba-split subband (Eg) of the 2DEG. This is 
entirely consistent with the results of spin-resolved pho- 
toemission for these statesP While it is difficult to re- 
solve the outer branch of the second Rashba-split sub- 
band (Ei), its spin-split pair is visible as the innermost 
Fermi surface sheet at positive k x values in Fig. S3 ^b-e) . 
The asymmetry of this state has the same sign as the 
inner branch of the lower subband state (Fig. S3 f)), in- 
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FIG. S3: Circular dichroism measurements of the Fermi surface, performed using the experimental geometry shown in (a). 
Fermi surface map, measured using (b) right-hand and (c) left-hand circularly polarized light, respectively. These give the 
dichroism and sum shown in (d) and (e), respectively, (f) Measured asymmetry (defined in the text) of the topological, 
Rashba-split inner and outer branch of the lower 2DEG and innermost branch of the upper 2DEG states. These are vertically 
offset for clarity, and fit by Asm(4>) + Bsin(30) (solid black lines), with the coefficients given in the figure. 



dicating an overall up-down-up-down-up ordering of the 
spin for consecutive Fermi surface pockets. 

The asymmetry of the inner branch of both the lower 
and upper (Eq and E%) subbands is well described by a 
single sinusoid, with negligible sin(30) component in the 
fits. The spin therefore lies almost entirely within the 
surface plane for these near-circular Fermi surfaces, as 
expected for a simple model Rashba systemPS On the 
other hand, the asymmetries of the outer branch of the 
lower (Eq) subband and the TSS both contain a substan- 
tial component which varies as sin(3(/>), indicating a sig- 
nificant canting of the spin out of the surface plane. This 
is entirely consistent with th e strong hexagonal warping 
of their Fermi surfaces! 12 * 15 * ^ 



At energies close to the band bottom of the 2DEG 
states, such an out-of-plane spin polarization is still 
present for the TSS (Fig. |S4[ c)). Now, however, the spin 
lies within the surface plane not only for the inner but 
also the outer branch of the lowest 2DEG subband. At 
higher binding energies still, well below the band bottom 
of the 2DEG states and approaching the Dirac point, the 
ideal in-plane helical spin texture is recovered also for the 
topological state (Fig. |S4|f)). 

Therefore, not only the direction of the spin winding, 
but also its effective dimensionality, has a distinct vari- 
ation with binding energy for each of the multiple elec- 
tronic bands that make up the complex surface electronic 
structure of these materials. 
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FIG. S4: Constant energy contours from 0.15 eV below the Fermi level performed using (a) right-hand and (b) left-hand 
circularly polarized light, respectively, (c) Extracted asymmetry and model fits for the topological and lower 2DEG states, 
vertically offset for clarity, (d-f) Equivalent measurements for 0.4 eV below the Fermi level, where only the topological state is 
present. 
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